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We discuss heavy quark production in deep inelastic scattering (DIS). 

Consider the effect of a charm quark (with mass to) being produced in 
neutral current DIS. For simplicity we concentrate on virtual photon exchange 
and neglect the bottom quark. 

Precise data are now available on F 2 (x, Q 2 ). In pQCD the operator prod- 
uct expansion and asymptotic freedom allow one to calculate the scaling vio- 
lations. This is incorporated in the factorization theorem where short distance 
effects are separated into coefficient functions and long distance effects into 
parton densities (for lowest twist terms), whose evolution is governed by the 
DGLAP equations. A mass factorization scale \i is thereby introduced. The 
corrections to such an analysis are of order A 2 /Q 2 . The coefficient and split- 
ting functions have been calculated to order a 2 with massless quarks. The 
first problem is how the presence of a massive c quark affects the above anal- 
ysis? In principle there should not be any problem because massless quarks 
were only used for calculational convenience. At small scales fi <C rn the three 
flavour splitting functions and parton densities should remain the same. The 
coefficient functions then contain terms in to. The c quark contribution in 
DIS F2,c,exact(x, Q 2 , m 2 ) can then be calculated in this fixed three-flavour 
number scheme with parton densities for u,d and s quarks (and antiquarks) 
and added to the light quark component i*2,LrGHT(#j Q 2 ) to form Fz(x, Q 2 ). 

As Q 2 increases large terms in ln(Q 2 /m 2 ) develop which should be re- 
summed to give stable predictions. This means that for Q 2 ^> m 2 the c 
quark becomes effectively massless and should be removed from the coefficient 
functions into the light parton densities. When this has been done there are 
four-flavour densities (including one for a c quark) and four-flavour massless co- 
efficient functions. i*2,ZM-VFNS [x, Q 2 ) is the zero-mass variable flavour number 
scheme representation for DIS. Since the evolution in the DGLAP equations 
covers all scales there should also be a variable flavour number scheme (VFNS) 
which interpolates between the two descriptions above. Clearly the schemes 
are complementary in that there are overlap regions where one scheme is as 
good as another and maybe much simpler to use, but one has to be explicit 
in the construction of a VFNS by specifying which parton densities and which 



coefficient functions are used. The ACOT scheme is an order a s VFNS and 
recently a .simplified version of it has been proposed a. The CTEQ5 parton 
density seta contains both fixed flavour densities and VFNS sets which use the 
ACOT scheme. 

What is known in higher order pQCD? Ini we established an all-order re- 
lation between I*2,ZM-VFNS an d ^2,c, exact which implied an all-order relation 
between light and heavy quark densities. The operator matrix elements con- 
taining the terms in powers of ln(/i 2 /m 2 ), which should be incorporated into 
boundary conditions on the evolution of the four-flavour densities (including a 
c quark density,) , were calculated to order a 2 . Starting from a three-flavour set 
of densities inQ we have constructed a set of four-flavour densities in El which 
satisfy these boundary conditions and resum the large logarithms. Each set 
has a different gluon density and satisfies the momentum sum rule. Thorne and 
Roberts Q tried to construct coefficient functions for use with such ZM-VFNS 
parton densities but concluded that the scheme was too complicated. They 
then proposed a different scheme based on continuity of dFz(x : Q 2 )/dln(Q 2 ) 
across flavour thresholds. This TR scheme has been implemented in order a s 
in the MRST98Q (and subsequent MRST) parton density sets. 

Two VFNS's now exist in order a 2 . One is called the CSN scheme i. It 
uses heavy quark coefficient functions rSwd begins with a term in order a®. 
The second is called the BMSN schemeBB and is based on massless coefficient 
functions. The differences between all the VFNS's can be attributed to three 
ingredients entering their construction. The first one is the mass factorization 
procedure carried out before the large logarithms can be resummed. The sec- 
ond one is the matching condition imposed on the c-quark density, which has 
to vanish in the threshold region of the production process. The third one is 
the construction of both three-flavour and four-flavour sets of parton densities, 
respecting the relations inB. Although designed to calculate F2(x,Q 2 ) both 
the BMSN and CSN schemes divide it into two pieces called F2. c (x, Q 2 , A) and 
F2,light(£, Q 2 , A) respectively. The quantity A separates the collinear c — c 
pair production from the noncollinear pair production and is the value of the in- 
variant mass above which they can be detected. This split makes F2. c (x, Q 2 , A) 
collinear safe. The BMSN and CSN schemes are designed so they yield the 
same result as the three- flavour NLO expression for F2 iC (x, Q 2 , A) at Q 2 = m 2 . 
In this region the four-flavour scheme fails. As Q 2 increases there are differ- 
ences between the BMSN, CSN and three-flavour schemes reflecting both how 
the higher order terms are resummed into densities and the influence of terms 
in m 2 /Q 2 . At large Q 2 both the BMSN and CSN results for F 2 , c (x,Q 2 ,A) 
agree numerically with the four-flavour results. For x » 10~ 3 and Q 2 = 10 3 
(GeV/c) 2 both disagree by approximately ten percent with the three-flavour 
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result showing that the resummation of the large logarithms is not a dramatic 
effect once one gets to order a 2 . The differences increase at larger x which, for 
fixed Q 2 , is closer to the threshold at Q 2 (x^ 1 — 1) = 4m 2 because effects due 
to terms in to are enhanced. 

Note that we have designed the soft/collinear split above without intro- 
ducing a fragmentation function. We can do this for suitably defined inclusive 
quantities. The next problem is the calculation of the single particle inclusive 
rate for producing D* (D*) mesons. Every experimental measurement must 
specify thresholds (or regions in phase space) below which the mesons cannot 
be detected. Under these circumstances potential collinear divergences remain 
in the final state. Is it better to keep the c-quark massive to avoid them or 
to take the c-quark massless and absorb the singularities into fragmentation 
functions (and resum them via evolution)? Either way we need additional 
information on the parametrization of the fragmentation functions describing 
the transition from quarks to mesons. Also there is another scale pt and a 
new set of logarithms. Which logarithms are more important? There is no 
gain in resumming terms in ln(Q 2 /m 2 ) into either the parton densities or the 
fragmentation functions if they are not dominant. Does one need a VFNS for 
single particle inclusive production or only a fixed flavour scheme and where 
do they provide a better description of the data? 

To begin we assume the differential production rate is measured at small 
Q 2 and is_modellcd by extrinsic production. The NLO exclusive program 
HVQDISEIJ uses three-flavour parton densities in the initial state and a massive 
c-quark to predict the rate. The reason for this is that potential collinear 
singularities involving the final state c-quarks are now controlled by m. We 
assumed a simple Peterson et al., type fragmentation function to model the 
transition from the c(c)-quark to the D* (D*) meson. The final parameters are 
the pole mass to, the fragmentation parameter e, the scale ^, which appears in 
the running coupling (and the parton densities) and the probabilitiz_for c — > D*. 
e differential rates for the D* (D*) mesons measured in the HlM and ZEUS 
experiments are in good, agreement with this NLO calculation. A recent 
comparison is available inEj. The experimental groups have used HVQDIS to 
integrate over all phase space to produce a quantity called F2 iC c(x,Q 2 ), the 
charm component of F 2 (x,Q 2 ). However note this is based on the HVQDIS 
program, which does not include diffractive production, intrinsic charm or 
resummation effects at large and/or small x. It is also not the same quantity 
we defined above with the parameter A. Authors who claim agreement between 
their theory and this experimental result should be aware of this fact. HVQDIS 
has been used to produce distributions in the invariant mass of the c — c (or 
D* — D*) pair and experimental results should be available soon. 



3 



The NLO calculation is complicated. Can one fit the same differential dis- 
tributions and/or the experimental i r 2,ce(^ J Q 2 ) in a- four-flavour scheme with 
a massless c-quark? This can be tried at two levels. First can one calculate 
the single particle inclusive distributions with massless quarks? Most of the 
D* (D*) data is at small pt so the only large scale is Q 2 (and fi 2 ). Therefore one 
needs to understand how to partition the potentially large terms in ln(Q 2 //i 2 ) 
and ln(/x 2 /m 2 ) between the parton densities and fragmentation functions leav- 
ing a collinear safe partonic differential coefficient function. Then how about 
F2,cc(x, Q 2 )? Can one integrate over the final states to remove the fragmenta- 
tion function and give a four-flavour description at the level of the coefficient 
functions? The latter investigation is more involved because every time we 
integrate over a variable we encounter a potential new collinear singularity as 
ttl — ► 0. These questions remain to be investigated. 

We end with the comment that the all-orders-, proof of factorization in 
DIS including heavy quark effects given by Collins E-3 contains a unitarity sum 
over all the final states and he therefore never defines a collinear safe charm 
contribution. Also the advantage of using massless coefficient functions in 
VFNS's was first proposed inH. 
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